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We determine transition probabilities in two exactly solvable multistate Landau-Zener (LZ) models and dis- 
cuss applications of our results to the theory of dynamic passage through a phase transition in the dissipationless 
quantum mechanical regime. In particular, we show that statistics of particles in a new phase demonstrate scal- 
ing behavior. Our results also reveal a symmetry that we claim is a property of a large class of multistate 
LZ models, whose explicit solutions are not presently known. We support our arguments by direct numerical 
simulations. 



I. INTRODUCTION 

The multistate Landau-Zener (LZ) problem is to deter- 
mine transition amplitudes among N discrete states after time 
evolution from — oo to +oo in systems described by the 
Schodinger's equation with linearly time-dependent coeffi- 
cients |[l]: 



(1) 



where A and B are constant N x N matrices. Any such 
a system can be transformed by a time-independent change 
of the basis into its canonical form, in which the matrix B 
is diagonal and the off-diagonal matrix elements of A are 
nonzero only between eigenstates of B with different eigen- 
values. Eigenstates of the matrix B are called the diabatic 
states. They transfer into eigenstates of the Hamiltonian when 
time is approaching ±oo. Therefore, for a finite number of 
coupled nondegenerate diabatic states, the scattering matrix 
can be defined. Arbitrary element 5„„' of such an x 
matrix S is the amplitude of the diabatic state n' at i — > +oo, 
given that at t — > — oo the system was at the state n. The 
related matrix P, Pn^n' ~ IS'nn'P, is called the matrix of 
transition probabilities. 

Studies of multistate Landau-Zener models (MLZMs) were 
pioneered by the article of Majorana ^ in which, in addi- 
tion to an independent discovery of the Landau-Zener formula 
12 , 3 1, he showed that any solution of a spin- 1/2 problem in 
time-dependent magnetic fields, including the LZ-model, can 
be generalized to the solution for dynamics of an arbitrary 
spin in the same magnetic field. Hence, the Majorana's ar- 
ticle |2 1 provided the transition amplitudes in the MLZM with 
the Hamiltonian H ~ (3tSz + gSx, where S is the arbitrary 
size spin operator Perhaps, the most unusual findings about 
the system ([T]i are the so called Brundobler-Elser formula [1 J 
and the no-go theorem |4| that provide expressions for some 
elements of the transition probability matrix in any model of 
the type ([T]i, suggesting that a strong progress in understand- 
ing complex MLZMs can be achieved. 

Today, there is no general recipe to determine all transi- 
tion amplitudes in a complex MLZM analytically. In order 
to achieve this, one would have to consider higher than 2nd 
order systems of differential equations with time-dependent 
coefficients. Nevertheless, a number of exact results lilllHSl, 



fully solvable nontrivial systems f^'-T^, and methods to study 
MLZMs 1 19. ,20il of the type ([TJ have been discovered. Some 
of the solved models were used in applications to condensed 
matter systems flT] . 

Here we add two models to the class of the worked-out 
solvable MLZMs. Our models correspond to transitions on 
a semi-infinite chain of sites with a possibility of pair-wise 
jumps between neighboring sites. Specifically, we will pro- 
vide transition probabilities for two models of type ([T]), whose 
Schrodinger's equations for amplitudes read: 

Model- 1: Square-Root Growing Coupling 



ia„ = — /3nia„+7\Ar+Ta„+i+7-v/na„„i, neN. (2) 
ModeI-2: Linearly Growing Coupling 

id„ = /3nto„ + "f{n + l)a„+i + 7na„_i, n e N, (3) 

with constant parameters f3 and 7, and where A^ is the set of 
nonnegative integer numbers. 

The reason why both models are solvable is because their 
Hamiltonians are quadratic when they are written in terms of 
bosonic creation/annihilation operators, a and a+. Thus, if 
we identify states \n) with eigenstates of the boson number 
operator a+a, then the ModeI-1 is the Schrodinger's equation 
for state amplitudes with the Hamiltonian 



Hi = -I3ta+a + 7(0+ + a), 



(4) 



and ModeI-2 is reproduced from the evolution of two coupled 
oscillators: 



H2 = I3td+a + 7(a+6+ + db), 



(5) 



under the condition that the initial population of the mode d is 
equal to the initial population of the mode b. 

Both models, Q and (|5]l, have an important practical re- 
alization. They describe a dynamic passage through a phase 
transition in a system of a molecular Rose condensate interact- 
ing with cold atoms near the Feshbach resonance. The Hamil- 
tonian of this system is usually written as Il22l424l : 



H 



E 



'/3tb+b + 



1_ 
N 



h.c. 



(6) 
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where b is the diatomic molecule annihilation operator and 
cip cr is the annihilation operator of a single atom with mo- 
mentum p and spin cr; N is the number of atoms and 7 de- 
scribes the strength of the conversion of atoms into molecules 
near the resonance. Model- 1 appears in the limit when atoms 
are in a macroscopic condensate state so that their operators 
can be treated as constant c-numbers. It describes creation 
of a molecular condensate beyond the mean-field assumption 
for the molecular field b. In particular, it can be applied to 
the experimentally important case with zero initial number of 
molecules. Model-2 corresponds to the opposite process of a 
decay of a molecular condensate into atomic modes. In this 
case, one assumes that the system has initially macroscopic 
number of Bose condensed molecules that pass through the 
Feshbach resonance and split into pairs of atoms with initially 
close to zero populations. 

Both models Q and (|5]l have been studied previously for 
application to the transition through the Feshbach resonance 
124-26 1, however, the focus was either on the average num- 
ber of molecule s/atoms converted during the process or on the 
evolution from the ground state. In this work we will explore 
the exact statistics of the number of defects (i.e particles in a 
new phase) created during the evolution starting from an arbi- 
trary initial state, and stressing universality and scaling in the 
full probability distribution of the possible outcomes of the 
sweeping through a phase transition process. 

The structure of our article is as follows. In Sections II and 
III, we derive the state-to-state transition probabilities in, re- 
spectively. Model- 1 and Model-2. In Section IV, we discuss 
possible applications of our results in the theory of dynamic 
quantum phase transitions. Appendix A is devoted to the nu- 
merical study of the symmetry of the transition probability 
matrix that we initially observed in solutions of our models 
and then claimed that it is actually the property of all MLZMs 
in linear chains. We summarize our results in Conclusion. 



II. TRANSITION PROBABILITIES IN MODEL-1 

The Hamiltonian of the Model- 1 can be easily transformed 
into the Hamiltonian of a quantum harmonic oscillator with 
a time-dependent force. Such systems were very well stud- 
ied previously for various applications. For example, one can 
engineer the Hamitonian of such a system in an array of opti- 
cal couplers and, literarily, observe the evolution of transition 
amplitudes with "time" ||27]| . 

To simplify our notation, first, we reduce the number of 
parameters by rescaling: t — >■ t/y/]3, g = 7/%/?. Next, we re- 
move strongly oscillating phase factors att ±00 by chang- 
ing the basis 



(7) 



which does not change the transition probabilities. After these 
transformations, the Hamiltonian Q is 



(8) 



We will search for the solution of the corresponding 
Schrodinger's equation in the form of a coherent state ansatz: 



„0(t)+Q(t)a+ 



|0). 



(9) 



Substituting this vector in equation i'i = Hi^ and collecting 
separately c-terms near and a^'i>{t) we find a pair of 
equations: 



i^ = ge^*^^^a, ia = ge 



(10) 



With initial conditions a{—oo) — q;_oo and 0(— 00) — 0, 
they have solutions 



a{t) = a-00 -ig dtie 



(t){t) = -iga^oo / rfiie'*i/2 
^ — 00 
ft pt 
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(11) 
(12) 

(13) 
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00 J ~oo 



dh / dt2e"i/^""=/^. (14) 



For the time-evolution from minus to plus infinity, this 

gives us a(-l-cxi) — a_oo — 5\/27re"^/'' and 0(+oo) = 
— g^TT + ga-ooV2Tre~^^^'^. Hence, if the initial state is a co- 
herent state, 5't_>._oo — \a-oo), the state after the transition 
becomes 



+ C30 



xe 



-gV2Tre" 



|0). 



(15) 
(16) 



A specific case of an initial state without defects: 5'_oo = 
|0) corresponds to a^oo = 0, and Eq. (I61 gives the wave 
function in the form 



+ 00 



(17) 



which is a coherent state with a 

vn+ \ r\\ X-^OO / n 



-g^/^^e^""/^. From 

e""^|0) = 'Y^^^oi'^^) / ^f^^-W) ^ we can explicitly read the 
transition amplitudes between states and n as = 
g-s''r(_(^y/27r)"e™'^/^/\/nT, and transition probabilities are 
explicitly given by 



^-2-Kg' 



(18) 



The number of defects (i.e. molecules in the created molecular 



condensate) has the Poisson distribution ( 18 1 with the average 
number (n) = 2715^. Returning to the parametrization (j4|i, we 
find (n) = 27r7^//3, i.e. the average number of defects scales 
as (3^^ with the sweeping rate. Two additional consequences 



follow from the knowledge of the full distribution ( 18 1. First, 
the Poisson distribution has all equal cumulants, so that this 
scaling is valid not only for the average number of defects 
but also for their variance, skewness and so on. The second 
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observation is that the probabihty of not creating any defect 
depends exponentially on the inverse sweeping rate /3: 



Pq^o 



(19) 



Later, we will return to this observation and claim that this be- 
havior is, in fact, more general than the scaling of the number 
of defects with the sweeping rate through a phase transition. 

In order to determine transition probabilities between any 
pair of levels n and rt', we use the fact that any state \n) can 
be written as a superposition of coherent states |e'^) as 



\n) 



= 1/2 



2^ ./o 



-inS\J.e 



(20) 



Combining this with ( fTT] ) we find that the state vector that 
starts as I ri) transforms into 



n!e 



+00 



27r 



(21) 



X exp [g^p2^e-'''l^e''^ + (e^^ - g^J2■Ke"l^)a^'^ |0>. 

For the transition amplitude to the state with n! > n we 
find 



n'>n : Snn' = e^^^ (-\/27rge*''/*) 

n ,, 



X 



E 



— ' m\(n — m)\{n' — n + m)! 

m=o ^ ' ^ ' 



/4v''-«,/— (22) 
(-2V)' 



Respectively, for n' < nwe find 



n < n 



(-V2^.9)"'+"e"("'""^/'' 



X 



n'!(-27r.g 



2\— m 



-^-^ m!(n' — to)!(7i — m)! 



(23) 



Transition probabilities Pn^n' — jS'nn' P can be compactly 
written in terms of the known special functions: 



n > n 



P . / — 



3-2V„!(27rg2)"'-» 



n' - n : P„^„ = p-^-s^ {L„{27:g^))' , (24) 

n <n: P„^„' = —^^ \^L^, (27r.g ) j , 

where L„ is the n-th Laguerre polynomial and is the asso- 
ciated Laguerre polynomial. In Fig. [T] we compare predictions 
(24 1 with our direct numerical simulations of (|2]). Results ap- 
pear to be in a perfect agreement with each other. 

Interestingly, amplitudes of transitions Snn' and Sn'n are 
different only by a phase factor, which corresponds to the 
symmetry of the transition probability matrix: 



(25) 




FIG. 1. (Color online) The numerical check (blue dots) of the theo- 
retical prediction (the red curve connecting discrete points) based on 
Eq. p.4\ for (a) transition probabilities from the state to the state n 
at g — 0.55; (b) transition probabilities from the state 2 to the state 
n at g = 0.55. Time evolution is from t = —600 tot = 600. 



This fact is quite surprising considering that there seems to be 
no obvious symmetry in the evolution equation of the model 
that would lead to ( |25] l. In fact, our numerical simulations 
show that this property is generic i.e. it is found in all MLZMs 
in finite and infinite linear chains. We will discuss this prop- 
erty in more detail in Appendix. 

Some transition probabilities with lowest n and n' are ex- 
plicitly given by 



Pi^i 

Pl^2 
P2^2 

Pl^3 
P2^3 



P2^i = e-^-^\ng'{g\-ir, 
e-2-9'(l-4V + 27rV)', 
P3^l = e-2'^«'^v^V(3-2V)^ 



(26) 



1 



(3 - 187rg2 + IgTr^g" - ^Tr^g^Y' 



It is useful to compare this list of transition probabilities with 
the ones for other known MLZMs. Thus, in all known ex- 
act solutions with a finite number of states, e.g. |2, 9, 10|, 
transition probabilities can be expressed as finite polynomials 
of exponents exp(— 7rg|^ //?ij), where gij and f3ij are, respec- 
tively, off-diagonal couplings and slope differences between 
pairs of diabatic states i and j. In another extreme, the ex- 
act solution of the MLZM in an infinite chain with constant 
couplings ifTSl returns transition probabilities in terms of the 
Bessel function: Pn^n' = Jfn-n'\iV^^^)- Our Model- 1 
shows features of both these classes. Transition probabilities 
contain an exponent e ^ but, unlike the known solved fi- 
nite size MLZMs, this exponent is multiplied by polynomials 
of 7rg^ rather than exponents of this combination. 



m. TRANSITION PROBABILITIES IN MODEL-2 

In order to solve Model-2, we adopt a different strategy 
ifTSll . Let us introduce the amplitude generating function 
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u = u{z,t) = y]„ fln^:". After rescaling, t — t/y/P and 
9 = l/^fP^ Eq. (pjl l3ecomes a linear partial differential equa- 
tion in terms of u(z, t): 



dtu + i[zt + g{z'^ + IjJO^zi = —igzu, 



(27) 



which can be solved by the method of characteristics, leading 
to equations: 



du 



-igzu, 



dz 
~di 



i[zt + g(z^ 



1) 



(28) 



(29) 



Nonlinear Eq. ( 29 1 can be transformed into the linear 2nd 
order differential equation by a change of variables 



z{t) 



idtajt) 

gait) 



where 



ita' — g^a = 0. 



(31) 



One of the initial conditions in ( [3T] l can be chosen arbitrar- 
ily. We will assume that |a(— oo)| = 1. The second initial 
condition is given by 



z_oo = a'(-cx))/a(-oo), 



(32) 



where 2_oo is the value of z{t) at t — > — oo. Substituting (30 1 
into (|28ll we also find 



u{t) = u^aca{t). 



(33) 



Equation (31 1 has two solutions with leading asymptotic s 

at t — > — oo: 



ai(-oo) - f9 



a2(-oo) 



-zi-'f -le'* (34) 



The Wronskian of such two solutions is equal to unity: W = 
ai{t)a'2{t) - a2{t)a[{t) = 1. Substituting (l34| into (|32| we 



find that 



Hence, 



a{t) = ai{t) + z_oo(5/i)«2(0- 
a[{t) + z^ao{g/-i)a'2{t) i 



z{t) 



ai{t) + z^ao{g/i)a2{t) g' 



(35) 



(36) 



In particular, this allows us to express z_oo in terms of z+c 
that is the value of z{t) at +oo: 



_ ai(+oo)z+oo(g/z) - a[{+oo) i 



a^(+oo) - z+oo(ff/i)a2(+oo) g' 



Assuming that the system starts on the level \n), we have 
U-oo = ^-oo- Substituting this into ([33]l we find 

(38) 



u{t) = z1^[ai{t) + z^^{g/i)a2{t)]. 




FIG. 2. (Color online) The numerical check (blue dots) of the the- 
oretical prediction (red curves) based on Eq. l[4TJ for (a) transition 
probabilities from the state to the state n at g = 0.5; (b) transi- 
tion probabilities from the state 2 to the state n at g = 0.4. Time in 
numerical simulations changes from t = —700 to t = 700. 



(30) From ( 38 i and (pTji, we obtain for z = z+ 



u{z) = 



ai{+Go)z{g/i) — a']^(+oo) i 
4(+oo) - z{g/i)a2{+oo) g 
1 



a2(+oo) - z{g / i)a2{+(x) ' 



(39) 



Asymptotics of ( [3T] i at t ^ +oo are known |23|. Their non- 
vanishing contributions read: 



ai(+oo) 

a2(+00) - f3'e~"''Sn^g')g-^^/4^/^;2^^g^:^/g^ 

4(+oo) - -ie'^f'e"'/^^-'^'. (40) 

It is possible to simplify the generating function u{z) by 
noticing that the multiplication of z by any complex number 
with a unit absolute value is equivalent to a mere phase trans- 
formation of all amplitudes a„, which does not influence the 
transition probabilities. Similarly, multiplication of u by an 
arbitrary phase factor, only shifts phases of a„ but not their 
absolute values. The time-dependent factors, e'* /^i^f , can 
be moved in front of u and the remaining time-dependent 
phase can be absorbed by z. Same can be applied to all factors 
g±jargr(jg-)^ 1, 1/i and e^^^^ in (40i. The remaining trimmed 



generating function of transition amplitudes from the state |n) 
to any other state then explicitly reads: 



m(z) 



z - Vl - 



-irg 



1 - zVl - / 1 - zVl - e^^Trg^ 



(41) 



In particular, by setting z = in ( [4T] ) and taking the square 
of the result we find the probability of a transition from the 
state |n) into the empty state |0): 



(42) 



which is the geometric distribution. From the previous studies 
of this model |26| we know that each defect that was present 
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at the initial state creates, on average, an exponentially large 
number of new defects, (n(+oo)) ~ n{—oo)e'^'^^ . Equation 



(42 1 shows that despite this tendency the transition probability 
into the empty state decays relatively slowly with the initial 
number of defects. Also, the result appears to be similar to the 
known solutions of MLZMs with a finite number of states in 
the sense that it depends only on the simple powers of e . 



According to (42i, the probability of creating no defects if 
the evolution starts at the ground state is given by a simple 
exponent: 



(43) 



where we again recalled that g — ^/ \fli in order to highlight 
the exponential dependence of Pqo on the inverse sweeping 
rate which is usually a directly controlled parameter in 
experiments. 

If the evolution starts with the ground state |0) then 



1 



1 



(44) 



Expanding (44i as a Taylor series we obtain the individual 
transition amplitudes. Taking their absolute square values we 
obtain transition probabilities from to all states in a closed 
form: 



(45) 



This particular distribution ( 45 1, but not the full amplitude 



generating function ( [4T] i, was previously found in |24| by a 
different approach. Comparing (45i and (42i we again find 
the symmetry P„^o = P^^n- 

Figure |2 shows a perfect agreement of the analytical pre- 
diction of (4T]i with our numerical simulations. Below, we 
also provide some explicit transition probabilities with lowest 
initial and final state indexes: 



Pl- 


1 = 




2e-2-9' 






Pl- 


2 = 


P2^i = e- 






-3e-2-9'))2, 


P2- 


2 = 






+ ee-^'^s' 


(46) 


Pl- 


3 = 


^3^1 = e- 


2^9^(1 _ 




4e-W)2^ 


P2- 


3 = 


P3^2 = e- 


2^9^(1 _ 




+ 


+ 2e 


-27rs 










^3- 


3 = 




2e-2-9' 


)2(l-10e- 


-2Trg'^ _ g-27rg^^-j2^ 



Again, we always find the symmetry of the transition proba- 
bility matrix: Pn^n' = Pn'^n- 



IV. APPLICATION TO DYNAMIC PHASE TRANSITIONS 

Modern studies of non-adiabatic effects during a passage 
through a quantum phase transition often concentrate on the 
almost adiabatic regime, in which a system almost reaches 
the new ground state after the sweep of a control parameter 



through a phase transition f28l. Nonadiabatic corrections in 
such a process are substantial only for small energy excita- 
tions defined by the sweeping rate /3. For systems with a 
continuous spectrum they result in scaling of the number of 
excitations (defects) with /?. Such a nearly adiabatic regime is 
usually hard to achieve in macroscopic and even mesoscopic 
quantum systems. There are many experimentally reported 
counterexamples to the naive application of the Landau-Zener 
formula to mesoscopic/macroscopic systems. Thus, an adi- 
abatic sweep of a magnetic field through the paramagnetic 
phase of an ensemble of initially polarized dipole-coupled 
nano-magnets usually leads to demagnetization rather than the 
transition to the ground state with the opposite polarization 
ll29l . Similarly, the conversion efficiency of an atomic con- 
densate into molecules after the passage through the Feshbach 
resonance is often close to 1 /2 rather than 1 at the lowest tem- 
peratures in the adiabatic limit ll30l[3Tl . Such deviations from 
the Landau-Zener theory have been explained, in particular, 
by the sensitivity of the adiabatic limit of the Landau-Zener 
formula to various decoherence effects |32|, which are gen- 
erally present in mesoscopic and macroscopic systems, and 
by the breaking of the adiabatic approximation by many-body 
interactions Ii24ll33l . 

The models that we solved in previous sections correspond 
to a different regime of a relatively fast sweep throughout a 
new phase. In this regime, the initial phase is mostly pre- 
served during the process, however, the number of defects 
(molecules in Model- 1 and atomic pairs in Model-2) can be 
large so that collective effects in their dynamics play the role. 
Such a regime is much less vulnerable to decoherence |34| 
due to much smaller time allowed for decoherence effects 
to accumulate, as well as due to specific robustness of the 
Landau-Zener processes to noise effects in this limit |35 1. 

We found that the probability distributions of the number 
of defects in this regime can be distinctly nonclassical. Con- 
sider, for example, a classical stochastic model of Markovian 
independent particle creation with a constant rate. The num- 
ber of particles in such a process has the Poisson distribution 
with the mean (n) and an exponentially small probability of 
creating zero number of particles Pq^o ~ e~^"^ In con- 
trast, for Model-2, it has been shown in 1261 that if the pas- 
sage through the phase transition starts with no defects, then 
^ e^'^^ , while we showed here that Pq^o = e^^'^^ , 
i.e. Po^o ^ (^i)^' which is dramatically different from the 
stochastic particle creation case. For example, if on average 
we observe 1000 defects, then it would be sufficient to re- 
peat the experiment ^ 1000 times to observe an event without 



a defect creation, while it would take 



,1000 



trials to ob- 



serve such an event if particles were generated by a classical 
stochastic process without memory. Our exact result shows 
that this property of the enhanced transition probability to the 
empty state persists even when the initial state is already pop- 
ulated with a number of defects. 

In studies of the quantum mechanical regime of a fast pas- 
sage through a phase transition, the probability to create no 
defects should play a special role also for the reason of its 
universal scaling with the inverse sweeping rate Eqs. 
( 19l and (43 i show that the probability Po-i-o is the same in 
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two models with very different behavior of other transition 
probabilities. This fact is not a coincidence but rather a direct 
consequence of the Brundobler-Elser formula [1], which was 
initially found in numerical simulations 1 1 1 and which by now 
is a mathematically rigorously proved result in the multi-state 
LZ theory ||5] |6l. According to this formula, if the system 
starts its evolution from the ground state and evolves accord- 
ing to Eq. ([TJ in time from — oo to +oo, then the probability 
to remain in the same diabatic state is given by 



Po- 



exp — 27r 



N 

E 



no^\ 



(47) 



where 7oi is the coupling between the state |0) and the exci- 
tation state and where /3o and /?; are the slopes of diabatic 
levels of these states. The summation is over all diabatic mi- 
crostates of the system. If a new phase is passed by changing 
the control parameter (such as an external magnetic field act- 
ing on a system of spins) from — oo to +oo with some rate 
j3 and if diabatic energies (elements of the diagonal matrix B 
in ([TJ) of all diabatic states depend linearly on the control pa- 
rameter then we have /3o ^ A — ctif3 for any state i, where 
parameters do not depend on the sweeping rate. The sur- 
vival probability then satisfies a universal scaling law: 



InPo- 



(48) 



where D is a constant, D = Stt^^-^ \ joi\''^/ai. Measure- 
ments of this constant by measuring Po-s-o can reveal a useful 
information about the coupling of the initial ground state to 
its excitations at the phase transition. We would like to stress 



again that the scaling law ( 48 1 is expected to be truly universal 
and not restricted by the conditions of the theory of the adia- 
batic passage through a quantum critical point [28 1. It should 
be valid for any linear passage through a region with a phase 
transition either with or without an exact level crossing point. 
In mesoscopic systems, such as atomic Bose condensates, the 
probability Pq^q can be made sufficiently large for measure- 
ments by increasing the sweeping rate f3. 
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FIG. 3. (Color online) Numerically obtained transition probabilities 
between the diabatic state n = 1 and other diabatic states in a 4-state 
MLZM in a linear chain. Off-diagonal couplings are proportional 
to the parameter g but with different, arbitrarily chosen coefficients: 
gi2 = 921 = 0.32g, 323 = 532 = 0.555, 534 = 543 = 0.75. 
Sweeping rates are /3i-/32 = 1.1765, /32-;53 = 1.5385, /?3-/34 = 
0.704. Results confirm the symmetry hypothesis. 
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FIG. 4. (Color online) Numerically obtained transition probabilities 
between the diabatic state n — 2 and other diabatic states in a 4-state 
MLZM in a linear chain. Parameters are the same as in Fig. [3] 



V. CONCLUSIONS 

We determined state-to-state transition probabilities in two 
multistate Landau-Zener models, which have practical rele- 
vance for the theory of dynamic passage through a Feshbach 
resonance. We highlighted the importance of the probabil- 
ity to generate zero number of defects by showing that it can 
be exponentially enhanced in quantum systems in compari- 
son to classical memoryless stochastic processes and that it 
shows the universal scaling InPo^o ~ as the function of 
the sweeping rate /3 through the phase transition. Our solu- 
tions also revealed the symmetry of the transition probability 
matrix that we explored numerically in Appendix for a set of 
finite chain models. All our tests support the hypothesis that 
in models of Landau-Zener transitions in chains the transition 
probability matrix is symmetric. 
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Appendix A: Symmetry of transition probability matrix of 
Landau-Zener models in linear chains 

Models 1 and 2 belong to the class of MLZMs that can be 
defined as Landau-Zener transitions in chains. Generally, the 
evolution for amplitudes of states in a chain can be written 
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FIG. 5. (Color online) Numerically obtained transition probabilities 
between the diabatic state n — 3 and other diabatic states in a 5-state 
MLZM in a linear chain. Off-diagonal couplings are proportional to 
parameter g but with different, arbitrarily chosen coefficients: gi2 = 
521 = 0.32g, = 532 = 0.55g, 534 = 543 = 0.7g, 545 = 554 = 
O.6I5. Sweeping rates are /3„ = n,n £ (1,2, 3, 4, 5). 
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FIG. 6. (Color online) Numerically obtained transition probabilities 
between the diabatic state n — 3 and some other diabatic states in a 
6-state MLZM. Off-diagonal couplings are proportional to parameter 
g but with different, arbitrarily chosen coefficients: 512 = 521 = 
0-855, 323 = 532 = 0.555, 534 = 543 = 0.75, 545 = 554 = 0.925, 
556 = 565 = 1-05. Sweeping rates are /3i — /32 = 1-2, /32 — = 
0.65, ^3 - ^4 = 0.7, /34 - A = 0.91, ft - /Sa = 1.0. 



as 



lai 



(Al) 

/3itai -I- 3^202, iaN = PNtaN -f gM-i.Na-N^i, 



where n e (2, . . . , — 1) with some real constants /3fe, k E 
(1, . . . N) and complex constants Qk.k+i, k £ {1, . . . N — 1) . 

Our explicit solutions of Model- 1 and Model -2 revealed the 
symmetry of the transition probability matrix: 



P. 



(A2) 



which is valid for transitions between any pair of diabatic 
states n and n'. This observation is surprising because the 



FIG. 7. (Color online) Numerically obtained transition probabilities 
as functions of time in a 4-state MLZM. Blue ciin>e: time dependence 
of the probability of being at the level n = 1 if the evolution starts at 
the level n — 2. Time evolution is from t = —700 tot — 700. Red 
curve: time dependence of the probability of being at the level n = 2 
if the evolution starts at the level n = 1. Black dots: The probability 
of being at the level n — 2 during the time reversed evolution from 
t = 700 backward in time. The evolution starts at the level n — 1 
at t = 700. Parameters of the model: 512 — 0.82, 523 = 0.55, 
534 — 0.7, Pn — n, n £ (1, 2, 3, 4). The symmetry Pi^2 — -P2-+1 
is found only asymptotically but at intermediate time, P2^i(i) 7^ 
Pi^2{t)- Time reversed evolution shows that P2^i (t) 7^ Pi^2{—t) 
but Pi^2{t) = Pi^2{-t). 



original equations that define our models explicitly break the 
chiral symmetry. Moreover, from some of the exactly solv- 
able models, such as the Demkov-Osherov model [16|, it is 
known that generally there is no such a symmetry in the full 
system ([T]i. On the other hand, in addition to our models, the 
available solutions of finite size MLZMs in linear chains that 
include the Majorana's solution for an arbitrary spin |2] and 
the exact solution of the 3-state MLZM |13| also show the 
symmetry ( |A2| i. Hence, based on the available information, 
we suggest the "symmetry hypothesis" that Eq. ( A2 1 is actu- 
ally asymptotically (i.e. for the time evolution from —00 to 
+00) exact for all MLZMs in linear chains. 



Today, the explicit expressions for transition probabilities 
in a general model (Al 1 are unknown. Therefore, to test the 
symmetry hypothesis, we performed a number of numerical 
tests with models of different complexity. Figs. |3|4|5|6| show 
some of our results for 4-,5-, and 6-state models of the type 
(Al I. Numerical tests with several other systems, including 
semi-infinite chains, of the type (All were also performed but 



not shown here. We found that all our numeiical tests sup- 



ported the hypothesis ( A2 1 in MLZMs of the type (All. 

The symmetry (A2i is not a simple consequence of some 
trivial symmetry of (Al 1. An example of a trivial symmetry 



is the statement that in any system ( |A1[ ) transition probabil- 
ities depend only on absolute values of coupling constants, 
l3n,ri+i|. This follows from the fact that a time independent 
change of the basis a„ — > a„e"^" linearly changes phases of 
9n,n±i- It is always possible then to choose N phases so 
that all — 1 parameters gn.n+i become real. Another trivial 
fact is that transition probabilities depend only on differences 
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of the level slopes f3i — f3j, which follows from the fact that 
a uniform shift of all f3k by a constant (3, (3k ^ Pk + P for 
k e (1, • • • N), is merely equivalent to the time-dependent 
phase shift a„ — >■ anC^^'^ Such "trivial" symmetries man- 
ifest themselves not only asymptotically but also during the 



full time of the evolution of the system. In contrast, Eq. ( A2 1 
is generally not satisfied at intermediate times during the evo- 
lution, as we illustrate in Fig. [7] Generally, Eq. ( A2 1 is sat- 
isfied only in the limit t +oo when both transition prob- 
abilities saturate at the same value. Figure [7] also shows that 
in MLZMs in chains the time-reversed process produces the 
same probability matrix but, alone, this does not explain the 
asymptotic symmetry (|A2 



The symmetry of the probability matrix considerably re- 
duces the number of unknown parameters. The importance 
of models of the type (Al i for the theory of the Feshbach 



resonance has triggered research on developing approximate 
analytical and numerical approaches to solve such models 
ll24l [33l [36l . The exact symmetry ( |A2| l can be a useful tool to 
test such approximations and reduce the number of unknown 
functions during studies of statistics of defects. On the other 
hand, the existence of such a nontrivial symmetry suggests 
that there is a way to rigorously prove it. Such a proof should 
shed new light on the properties of MLZMs and, possibly, 
it may lead to the expression for the scattering matrix of the 



whole problem ( Al i in terms of the known special functions. 
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